New possible source of huge neutrino bursts 
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It is shown that primordial black holes (PBHs) of masses M > 5 x 10 14 g are able to absorb white 
dwarfs (WD) completely for the time of their existence, giving rise to formation of black holes of 
WD mases. The WD absorption is accompanied by up to W 52 erg neutrino bursts which could both 
considerably contribute to the cosmic neutrino flux and open up a new way of PBHs detection, 
especially effective in placing new constraints on abundance of nonevaporating PBHs with masses 
M > 10 15 g. 
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PBHs could have formed in the early universe ||] 
from initial inhomogeneities, topological defects (strings, 
bubbles) or at cosmological phase transitions. PBHs pro- 
vide a unique probe of the Big Bang and even their non- 
existence gives vital cosmological information |3|, ||, |i| . 

The Hawking radiation || allows to set up severe con- 
straints on the PBH abundance. In particular, the ra- 
diation of all the evaporated PBHs could have given its 
contribution to the diffuse cosmic flux [Q, |[. Very 
stringent bound has been obtained from the absence of 
a detectable diffuse flux of 100 MeV photons 0. It was 
claimed, however, that these bounds can be improved 
using available data on diffuse neutrino flux || |j| . 

Another way of PBHs search consists in attempts to 
detect gamma fjj or neutrino J§| bursts accompanying 
the final stage of PBH evaporation. The PBHs which 
complete their evaporation at present have had the ini- 
tial mass Mo equal to the critical one (Hawking mass) 
M* ~ 5 x 10 14 g. For all the time of their evaporation 
such PBHs could emit about M*c 2 ~ 5 x 10 35 er<?, where c 
is the speed of light, what is much less, for example, than 
the energy of supernova explosion reaching 3 x 10 53 erg. 
The comparatively low energy release undoubtedly con- 
strains the search possibilities of the most intensively 
evaporating PBHs and makes utterly ineffective that of 
weakly evaporating PBHs with M ~ Mq ^> M* , the con- 
straints on which are more than ten orders lower than 
that on PBHs with M ~ A/* |, @, §. 

We would like to suggest here an alternative method of 
PBHs search, free of this deficiency. Namely, we will show 
that being trapped inside a white dwarf (WD) a PBH of 
mass Mq > Af* can completely absorb it initiating a huge 
neutrino burst with the energy reaching 10 52 erg. 

First, we will consider the simplest case of nonrotat- 
ing cool WD and apply the general relativistic theory 
of spherical, steady-state, adiabatic accretion of contin- 
uous Pascal liquid [|l0| |ll| to show that the time of a 
complete WD absorption by a Schwarzschild PBH can 
be less than the age of the universe. A special considera- 
tion has shown that the heating caused by PBH Hawking 
radiation, WD matter compression and neutronization 
does not violate the high degeneracy of WD matter. The 



mass-energy density E and pressure P of the degenerate 
WD matter are given by the equations jl0|, |ll|, [l2| 
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where n is the electron number density, 
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is the electron Fermi momentum pf{p) = (37r 2 n) 1 / 3 fi. ex- 
pressed in m e c units, m e and Ti are the electron mass 
and the Planck constant, respectively, fj, is the WD mat- 
ter mass which falls at one electron. We will consider the 
most widespread case of carbon WDs for which /i ~ 2m„, 
where m n is the nucleon mass. The term of Eq.(l) con- 
taining the small parameter 
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where a ~ 1/137 and Z is the atomic number, describes 
the Coulomb correction flo|| which is principally impor- 
tant in the case of low WD density. 

The spherical, steady-state accretion flow of degener- 
ate Fermi gas is described by the radial component of 
4-velocity u and the electron number density n. Both 
of them depend only on the radial Schwarzschild coor- 
dinate measured from the PBH center and obey |n], [ll]] 
the electron conservation equation 



M = 47r/in(r)cw(r)r 2 = const 
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which gives the rest mass accretion rate as well as the 
relativistic Bernoulli equation 
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containing the gravitational radius r g — 2GM/c 2 ~ 
1.5(M/10 15 g)10~ 15 m where G is the Newton's constant. 
The steady-state accretion assumes that n — > no = const 
at r — > oo. The constant signs in Eqs. (5) and (6) reflect 
the independence on r. 

The basic point of the theory (loj [nj consists in ex- 
istence of a special (sonic) point r — r s at which the 
relations 

hold where a s = a(r s ) is the sound speed a = {dP/dE) 1 / 2 
expressed in c units and taken at r — r s . Using the Em- 
den solution (loj [l^] for the WD density distribution one 
can make sure that, if a PBH with mass M is situated at 
the WD center, the uniformity condition in the essential 
accretion region r ~ r s is fulfilled at M < 10 -3 Mwd 
where Mwd is the WD mass. The electron number den- 
sity no in Eq. (6) can be considered as that at the WD 
center in this case. 

Substituting Eqs. (1), (2) and (7) into Eq. (6) one 
obtains in the first order with respect to m e //i ~ 2.74 x 
1CP 4 the equation 
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connecting x s — x{n(r s )) with xq = x(no)- Eq. (8) 
demonstrates that the Coulomb correction is indeed es- 
sential at the density p = E/c 2 < 100g/cm 3 giving rise 
to the existence of two solutions of Eq. (8) at xq < 2a' 
(a' ~ 0.015 at Z = 6). However here we will consider 
the usual case of much denser WDs in which xq ^> 2a' 
and Eq. (8) has a single solution taking the simple form 
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or p{r s ) ~ 8po m the ultrarelativistic limit of xq 3> 1, 
po ^> 2 x 10 6 g/cm 3 . Eqs. (5) and (6) also allow one to 
calculate the parameter (3) maximum value 

(10) 



at r = r g . The corresponding Fermi energy Sp(r g ) ~ 
m e cx(r g ) reaches 20 — lOOMeV. 

Substituting x s into Eq. (7) one finds the values n{r s ), 
r s and u(r s ) allowing to calculate the right hand side of 
Eq. (5) and to obtain the law 
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of the PBH accretion mass growth starting from its initial 
value Mq = M(0) and characterized by the time 
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which can be considered as that of WD absorption. It 
should be emphasized that though Eq. (11) loses its ap- 
plicability at M(t) > 10~ 3 M\yd, the accuracy of the 
estimate (12) is quite high since Eq. (11) does not de- 
scribe only a very small part At/T a i, s — Mo / (10~ 3 M\y d) 
(10~ 16 -th part at M - M*) of the time (12). Thus, the 
approach [ 10[ [ll[] has led us to a very precise estimate of 
the time of WD absorption by a PBH. 

Using the relation (8) of x s with xq = x(no) one can 
calculate the time (12) dependence on the WD central 
density po = E(no)/c 2 . Fig. 1 demonstrates that a 
PBH with initial mass Mo = 10 15 <7 can absorb rela- 
tivistic (po > 10 6 g/cm 3 ) WDs for several Gyrs. Since 
T abs oc Mq 1 , PBHs with initial masses M > 10 17 g will 
absorb WDs of any density for the time of their exis- 
tence. The absorption process could give rise to forma- 
tion of black holes of WDs masses which can not form in 
the star collapse process. Since WDs are quite abundant 
their absorption by PBHs could, in principle, be more 
productive source of star-mass black holes than the star 
collapse. 
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FIG. 1: Time of a complete WD absorption by a PBH with 
mass Mo = 10 1 g as a function of WD central density. 

A WD absorption by a PBH would be observed as its 
sudden disappearance accompanied by a moderate light 
(x-ray) flash. Such a disappearance, however, is difficult 
to observe at distances exceeding lOOps. We will show 
now that a WD absorption by a PBH is accompanied by a 
huge neutrino burst detectable at much longer distances. 
Neutrinos will be produced in neutronization reactions 
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Be + v in carbon WDs 



in the region with electron Fermi energy £f(t) 



rrieC 2 y/1 + x 2 (r) exceeding the neutronization threshold 
energy A z (A c = 13.88MeU, A B = 12.17MeU). Since 
£p (r g ) exceeds A c considerably at any xq and po (see 
Eqs. (8) and (10)) a PBH will give rise to the neutron- 
ization process in any carbon WD. 

Both the local neutronization probability of a nucleus 
with atomic number Z flOtl 
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where 

f'de, 

f = s(e 2 ~ miosis ~A z f, 
and the accompanying neutrino emission intensity 

ln2 
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II = Il{e F (r)) = 



r z .i/(m e c 2 , A z ) 
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are determined by the local electron number density n(r), 
the threshold energy A z and the half-life T z _ x of the 
daughter nucleus (T B = 20.20ms, T Bo = 21.31ms). 

Neutronization begins at the " neutronization radius" 
r n > r g at which £p (n(r n )) = A c and p n — E(n(r n )) / c 2 
= 3.90 x 10 10 g/cm 3 Q. According to Eq. (9) one has 
r n = r s at the central WD density p ~ p n /8 ~ 4.90 x 
10 9 (//cm 3 . Since most of WDs have much smaller ones 
let us assume that r„ < r s . The WD matter motion in 
the neutronization region r g < r < r n is determined in 
this case rather by gravity than by pressure and is well 
described by the free fall equation 



u ( r ) = \frj~ r 



(16) 



written in Lagrangian coordinates. Eq. (16) allows both 
to introduce a proper time r(r) = f n dr/[cu(r)] , r(r n ) = 
0, of elementary volume motion from r n to r and to write 
the electron number conservation equation r 2 nu — const, 
which holds in neglect of electron absorption, in the form 
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Eqs. (8) and (10) obtained in the model of steady-state 
accretion jl(| [ll[] fail to describe the accretion rate sat- 
uration accompanied by the most intense neutrino emis- 
sion. The last comes from the region r ~ r max giving the 
main contribution to the integral determining the WD 
mass. According to the Emden solution for ultrarela- 
tivistic Fermi gas |T^| r max equals 0.237i?^_o, where 
Rwd is the WD radius. The WD matter with the ini- 
tial density p(r max ) — 0.312po experiences the attraction 
of the mass 0.37 Mwd- Though this mass is contracted 
in the process of absorption by a PBH it determines the 
constant value of r g in Eqs. (16) and (17) describing the 
final stage of accretion accompanied by neutronization. 
The nonuniformity of the region r ~ r max giving the 
main contribution to the neutrino emission will be taken 
into consideration by using p(r max ) instead of po in Eqs. 
(8) and (10) determining the x max value. 

Let us introduce the C, B and Be nuclei fractions nor- 
malized by the equality c c + c B + c Bc = 1 and satisfying 
the initial conditions c c (0) = 1, c B (0) = c Bo (0) = at the 
moment r(r„) = of passage through the neutronization 
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FIG. 2: Characteristics of neutrino emission process accom- 
panying a carbon WD absorption by a PBH: (a) average neu- 
trino number which falls at one nucleus, (b) average neutrino 
energy and (c) neutrino burst energy as a function of cen- 
tral WD density as predicted by the model described in the 
text (dashed line) and by a more sophisticated approach (full 
lines). 



radius. The first two fractions obey the equations 



dc c /d,T = -r c c c 
dc B /dr = r c c c - r B c B . 
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Eqs. (16) and (17) allow to pass to the variable r in Eqs. 
(13) and (18) and easily integrate the last yielding the 
analytical expressions for c c (r) and c B (r), as well as for 
both the average neutrino number 



< N v >= 2 — 2c c (r 5 ) - c B (r g ) 
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which falls at one carbon nucleus and the average neu- 
trino energy 
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The < N u > and < e v > dependencies on the central 
WD density p are given, respectively, on Figs. 2(a) and 
(b) along with the same dependencies evaluated using a 
more sophisticated model taking into account the elec- 
tron absorption, the pressure influence on accreting mat- 
ter motion and the nonequilibrium heating of the last 
which accompanies the neutronization process [ p^[ . The 
negligible difference in predictions of these two models 
justifies the undertaken consideration of the simpler one. 

Fig. 2(a) demonstrates the possibility of complete car- 
bon and subsequent boron neutronization. The average 
and maximum neutrino energies can reach 30 (see Fig. 
2(b)) and lOOMeV (see Eqs. (9) and (10)), respectively. 
The dependence of the energy 
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of the neutrino burst accompanying the total WD ab- 
sorption on the central WD density is given on Fig. 2(c) 
demonstrating that E v can exceed 10 52 erg. The neutrino 
burst duration will not considerably exceed the free fall 
time, since neutrino will travel through the WD matter 
nearly freely. The huge neutrino flux, nevertheless, will 
heat the WD matter up to T > 10 7 K and, possibly, ini- 
tiate the nuclear burning |Q in the most dense WDs. 
These processes will be considered elsewhere. 

Thus, we have demonstrated that the energy of a 
purely neutrino burst accompanying a WD absorption by 
a PBH is only a few orders less than that of a neutrino- 
antineutrino burst accompanying a supernova explosion 
[ p~5|| . If such a burst occurs at the typical distance of 
lOkps it will cause up to 50 events of ve scattering in the 
Super-Kamiokande detector and, thus, can be detected 
with confidence. In fact, it looks feasible to design a 
detector able to detect most of the predicted neutrino 



bursts in our Galaxy g, A possible frequency of 

these bursts depends on such uncertain factors as PBH 
abundance and mass distribution as well as on the prob- 
ability of PBHs presence in WDs or capture by them. 
Taking into consideration the high WDs abundance in 
the Galaxy and assuming more or less uniform distri- 
bution of the moments of WDs complete absorption in 
time, optimistically one can hypothesize that the pre- 
dicted neutrino bursts can be detected nearly every year. 

According to Ref. §] PBHs with masses M < Ah 
contribute mostly to the diffuse neutrino flux, the study 
of which already allows to set up restrictions on the power 
of the "blue" spectrum of initial density fluctuations. 
The possibility of intensive neutrino emission by WDs 
being absorbed by PBHs radically change this picture. 
First, since the energy (21) 11-16 orders exceeds the rest 
energy M*c 2 , even a small number of PBHs absorbing 
most weakly rotating WDs can improve the bounds on 
PBHs absence by many orders. Second, since the heav- 
ier PBHs absorb WDs more readily, the diffuse neutrino 
flux, in fact, allows to set up the most stringent bounds 
on the nonevaporating PBHs with masses Mq M*. 

Thus, we have demonstrated that PBHs can absorb 
WDs for the time of their existence giving rise to forma- 
tion of black holes of WD's masses. The WDs absorption 
is accompanied by the neutrino bursts which can be de- 
tected directly and could have given a considerable con- 
tribution to the diffuse neutrino flux. Most considerably 
these neutrino bursts widen the possibilities of search of 
nonevaporating PBHs with masses Mq ^> 10 15 g. 
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